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ABSTRACT 

We study the observability of the reionization epoch through the 21 cm hyperfine 
transition of neutral hydrogen. We use a high-resolution cosmological simulation (in- 
cluding hydrodynamics) together with a fast radiative transfer algorithm to compute 
the evolution of 21 cm emission from the intergalactic medium (IGM) in several dif- 
ferent models of reionization. We show that the mean brightness temperature of the 
IGM drops from STf, ~ 25 mK to 8Tb ~ 1CP 2 mK during overlap (over a frequency 
interval Av ~ 25 MHz), while the root mean square fluctuations on small scales drop 
abruptly from (ST 2 ) 1 / 2 ~ 10 mK to {ST 2 ) 1 / 2 ~ 1CT 1 mK at the end of overlap. We 
show that 21 cm observations can efficiently discriminate models with a single early 
reionization epoch from models with two distinct reionization episodes. 

Key words: cosmology: theory - intergalactic medium - diffuse radiation 



1 INTRODUCTION 

One of the major challenges in cosmology is to under- 
stand the first phases of structure formation. We do not 
yet know whether early generations of sources were simi- 
lar to galaxies and quasars in the local universe or quali- 
tatively different from them (and if so, in what ways). Be- 
cause these sources affect the formation of all subsequent 
generations of luminous objects, understanding their char- 
acteristics is critical to any theory of structure formation. 
In many ways, the defining event of this era is "reioniza- 
tion," when the first sources of light ionized hydrogen in 
the intergalactic medium (IGM). The timing of this event, 
and its morphological evolution, contains a wealth of infor- 
mation about these first sources and about the IGM itself 
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Observational constraints on this epoch are difficult to 
obtain, but those that do exist are intriguing. The most 
straightforward way to study the ionization state of the IGM 
is through quasar absorption spectra: regions with relatively 
large HI densities appear as absorption troughs in the quasar 
spectra, creating a "Lya forest" of absorption lines. Spectra 
of z ~ 6 quasars selected from the Sloan Digital Sky Survey 1 
(SDSS) not only show at least one extended region of zero 
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transmission feecker et all feoOl) but also impl y that the 
ioniz ing background is rising rapidly at this time (iFan et alJ 
2002). Both of these results are consistent with a scenario 
in which reionization ends at z r ~ 6. Another constraint 
comes from measuring the fraction of the cosmic microwave 
background (CMB) radiation that has been Thomson scat- 
tered by ionized gas in the IGM. Recent observations made 
by the Wilkinson Microwave Anisotropy Prob e 2 ( WMAP) 
imply that reionization occurred at z r > 14 llKoeut et alJ 
l2003tlSpergel et alJ2003tl . Finally, measurements of the tem- 
perature of the Lya forest at z ~ 2-4 sug gest an order 
unity change in the ionized fraction at z r < 10 llTheuns et alJ 
2002; Hui & Haiman 2003), although this argument is sub- 
ject to important unc ertainties about Hell reionization (e.g., 
ISokasian etafl 12002). Reconciling these different observa- 
tions implies that sources must exhibit substantial (perhaps 
qualit ative) evolution between high redshifts and the present 
dav JSokasian et alj|2003t IWvithe fc Loebl2003bl: ICenll2003l: 
lHaiman fc Holderl2003l) ~A better understanding of reioniza- 
tion has the potential to teach us the detailed evolutionary 
history of luminous sources. 

Unfortunately, learning much more from these obser- 
vations promises to be difficult. The optical depth of the 
IGM to Lya absorption is r a « 6.45 x 10 5 :ehi[(1 + z)/10]' i/2 
jGunn fc Petersonlll965l) . where xm is the neutral fraction 
and where we have assumed the currently favored cosmolog- 
ical parameters (see below). A neutral fraction im > 10~ 3 

2 See http://map.gsfc.nasa.gov/ 
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will therefore render the absorption trough completely black; 
quasar absorption spectra can clearly probe only the very 
late stages of reionization. CMB measurements, on the other 
hand, depend on the ionized gas density and thus are most 
sensitive to the early stages of reionization. However, ex- 
isting observations essentially provide only an integral con- 
straint on the ionized gas column; distinguishing different 
reionization histories wit h CMB observation s alone promises 
to be quite challenging l|Holder et, aljEoO^I . 

It is therefore crucial to develop other probes of the 
high-redshift IGM. One possibility is to observe a transi- 
tion much weaker then Lya. Perhaps the most interesting 
candidate is the 21 cm hyperfi ne line of neutral hydrogen 
in the IGM (lFieldlll95SL Il959alb The physics of the 21 cm 
trans ition has been well-s t udied in the cosmological context 
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Miralda-Escudc 2003). There are three meth- 
ods to study the IGM through this line. First, we can 
search for emission or absorption from the diffuse IGM and 
condensed minihaloes (i.e., collapsed objects that are un- 
able to cool and form stars) at high redshifts. While the 
absolute signal from this gas will be swamped by fore- 
ground emission, the signal on the sky will fluctuate be- 
cause i^f^^iations^n th^^G^^ensity and neutral frac - 
tion jMadau et alll997tlTozzi et al.ll2000l : llliev et all2002fr . 
Second, we can seek a global (all sky) signature when the 
neutral emitting gas is destroyed. Assuming that reioniza- 
tion occurs rapidly, this will appear as a "step" in the 
otherwise smooth l ow-frequency radio spectrum of the sky 
JShaver et al]ll999h . Third, high-resolution spectra of pow- 
erful radio sources at high-redshifts will reveal absorption 
features caus ed by sheets, filaments, and minihaloes in th e 
neutral IGM dCarilli et alJl2002t iFurlanetto fc Loebll2002t) . 

To date, predictions about the 21 cm signal of the dif- 
fuse IGM near the epoch of reionization have been difficult 
to obtain. This is because such predictions require a care- 
ful treatment of ionizing sources and radiative transfer in 
order to describe adequately the transition from a neutral 
to an ionized medium. Only recently has it become possible 
to combine numerical radiative transfer schemes with cos- 
mological simulati ons in order to study this process in de- 
tail <GnediJ200fllRazoumov et alJl2002HCiardi et alJl2003t 
ISokasian et all2003ft . As shown bv lCiardi fc Madaul fcool) . 
such simulations allow us, for the first time, to realistically 
model the behavior of 21 cm IGM emission . In th is paper, 
we analyze the simulations of lSokasian et all <l2003f) in order 
to quantify the 21 cm emission expected before, during, and 
after reionization. We focus purely on the emission signal be- 
cause absorption studies require careful treatment of both 
ultraviolet radi ation (above and below the Lym an limit) and 
X-ray photons llChen fc Miralda-Escudll2003l) . which is ex- 
tremely difficult (see lCariUi^t^lTl2002T f'or a first attempt). 
We study scenarios in which the IGM is reionized once (at 
either early or late times) as well as "double reionization" 
histor ies that may reconcile the SPSS an d WMAP observa- 
tions JWvithe fc Loebll2003bl : ICenll2003h . We show that 21 
cm measurements offer a promising route for distinguishing 
reionization scenarios that would appear identical in CMB 
and quasar observations. 

We briefly describe our cosmological simulations and 
reionization models in &\ We then discuss the physics of 21 
cm emission, and our analysis procedure, in J|l We present 



our results in and discuss their observability and their 
implications for understanding reionization scenarios in jJSj 



2 SIMULATIONS 

All of our calculations are performed within the 5 cos- 
mological simulation of ISpringel fc Hernauisd i2003 b*l . This 
smoothed-particle hydrodynamics (SPH) si mulation uses 
a modified version of the GADGET code (ISpringel et alJ 
l200ll) incorporating a new conservative formulation of 
SPH with the specific entro py as an independent variable 
JSpringel fc Hernauisdl2002l) . The simulation also includes 
a new description of star for mation and feedback in the in - 
terstellar medium of galaxies (ISpringel fc Hernauisdl2003aD . 
Within the extended series of simulations performed by 
ISpringel fc HernauistH2003rJ) . this model yields a converged 
prediction for the cosm ic star formation rate . Using sim- 
ple physical arguments, iHernauist fc Springell J2003) have 
shown that the star formation rate density evolves accord- 
ing to: 
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Here the constants a — 0.012, = 0.041, and p*(0) = 
0.013 Mp, yr -1 Mpc~ 3 are chosen by fit ting the simulation 
results of ISpringel fc Hernauisd J2003bl) . 3 The form of the 
function is motivated by two competing processes affecting 
the total star formation rate. At high redshifts, when cool- 
ing times are short, star formation is only limited by the 
growth of massive haloes. At low redshifts, cooling times 
become long because of the expansion of the universe; thus 
in this regime the star formation rate density evolves pri- 
marily because of the expansion of the universe. 

The simulation assumes a ACDM cosmology with fl m — 
0.3, Ov = 0.7, n b = 0.04, H = lOO/i km s" 1 Mpc -1 (with 
h — 0.7), and a scale-invariant primordial power spectrum 
with index n = 1 normalized to as = 0.9 at the present day. 
These parameters are consis tent with the most recent cos- 
mological observations (e.g.. ISpergel et al1 l2003). The par- 
ticular simulation we choose has 324 3 dark matter particles 
and 324 3 SPH particles in a box with sides of Wh comov- 
ing Mpc, yielding particle masses of 2.12 x 10 6 /i _1 M Q and 
3.26 x 10 ft -1 Mq for the dark matter and gas components, 
respectively. 

The rad ia tive transfer is fully described in 
ISokasian et alJ i200ll) . so we outline only its most sig- 
nificant aspects here. We perform the adiative transfer 



3 We pause here to not e an error in Figure 12 of 
ISpringel fc Hernauisd i2003bl) . where the simulated star forma- 
tion rates were compared to observations. The observational re- 
sults shown in the Figure for < z < 2 should be multiplied by a 
factor h for comparison to the simulation. When this correction is 
applied, the observed star formation rate densities decrease and 
come into closer agreement with the simulation results; thus the 
simulation agrees with observations for all z < 2 about as well as 
is shown by the points at z = in the Figure. 
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on a Cartesian gri d with 200 3 cells usi ng an adaptive 
ray-casting scheme (I Abel fc Wandeltj 12002^ . Density fields, 
clumping factors, and source characteristics are taken 
from outputs of the hydrodynamic simulations. Thus, the 
radiative transfer is done in a post-processing step after 
the simulations have been completed, and we neglect the 
dynamical feedback of the radiation on structure formation. 
We turn off periodic boundary conditions for the radiation 
once the ionized volume exceeds a certain level. After this 
time we add photons that reach the edge of the box to a 
diffuse background. The calculation follows the reionization 
of both hydrogen and helium. Because ionizing radiation 
from small galaxies has a substantial effect on reionization, 
we select sources aggressively (including all groups of 16 or 
more dark matter particles). At each redshift, we correct 
the star formation rat es in low-mass haloes to ma tch the 
converged prediction of Hernquist fc Springe! (|2003). 

Because the source and IGM properties are determined 
self-consistently by the simulation, the major unknown pa- 
rameter in simulations of reionization is the efficiency with 
which each galaxy produces ionizing photons. This parame- 
ter includes two physical effects: the intrinsic rate at which 
newly formed stars create ionizing photons and the fraction 
of such photons that are able to escape the galaxy. As a 
fiducial value, we assume an intrinsic ionizing photon pro- 
duction rate of Ar i]fid = 10 53 (Af*/ M yr" 1 ) s" 1 , where M* 
is the star formation rate, corresponding to a Salpeter ini- 
tial mass function (IMF) with solar metallicity. We then 
set the actual production rate Ni = /ionA^fid- Any value 
/ion < 1 can be accomodated by the fiducial IMF, while 
/ion > 1 will require a change in the IMF. For example, 
very massive metal-free stars are ~ 30 more e fficient at cre- 
ating ionizing photons than "normal" stars dBromm et alJ 
l200lf ). ISokasian et alJ d2005t) found that the choice f ioa = 0.2 
matches data from z ~ 6 quasars selected th rough the SDSS 
jBecker et al J l200lt iDiorgovski et all 1200 lfl . However, this 
model implies an optical depth to Thomson scattering signif- 
icantly below the limi ts from recent observations by WMAP. 
ISokasian et al.l <|2003T) showed that one way to increase the 
optical depth without destroying the match to the quasar 
observations is to gradually increase the ionizing efficiency 
with redshift in the form f ioa (z) = 0.20e (z_6)/3 . Such a 
model achieves overlap at z ~ 13. However, because /i on 
declines rapidly with cosmic time, the statistics of quasar 
absorbers remain essentially unaffected (i.e., most of the ex- 
tra high-redshift photons are balanced by recombinations 
before z ~ 6). We note that in this model /i on > 1 for 
z > 11, indicating that the IMF at these redshifts differs 
from that of local galaxies. 

Another way to reconcile the quasar and CMB obs erva- 
tions i s either partial or compl ete "double reionization" iCenl 
2002: IWvithe fc Loebl l2003ai. In such a scenario, sources 
reionize the universe at z > 15. Some form of feedback then 
prevents these sources from continuing to form and grow, 
rapidly reducing the production rate of ionizing photons. 
For example, feedback can either prevent stars from form- 
ing in small haloes or change the mode of star formation. 
(In our notation, /i on essentially changes discontinuously at 
the moment of overlap.) Because the recombination time at 
high redshifts is short compared to the Hubble time, much 
of the universe recombines before "normal" protogalaxies 
can cause a second reionization event at a later time. Note 
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Figure 1. The integrated optical depth to electron scattering 
T e between the present day and redshift z. The curves corre- 
spond to our different reionization models. The solid curve as- 
sumes a constant /; on = 0.2 for all sources throughout cosmic 
time. The dot-dashed curve assumes a variable fi on (z) (see text). 
The short-dashed and long-dashed curves assume /i on = 0.2 but 
set xhi = 0.01 everywhere until zt = 14.5 and zt = 18.4, respec- 
tively. The horizon tal dotted lin e shows the lower limit to r e (at 
68% confidence) of lKogut et al] <2003ft . 



that, although both smooth and discontinuous increases in 
/ion produce large optical depths to the CMB, the qualita- 
tive evolution of the neutral gas differs markedly between 
the two models. Unfortunately, our simulation cannot re- 
solve the haloes that are responsible for the first reionization 
epoch in these scenarios. We therefore introduce an approx- 
imate model for double reionization by setting the neutral 
fraction xm = 0.01 throughout the simulation volume for 
all z > Zt' At z = Zt, we allow the gas to recombine. We 
set /ion = 0.2 for these two models. Note that the simula- 
tion assumes that the ionized gas has T = 5000 K. If the 
photoionized gas has a higher initial temperature, recom- 
bination may take slightly longer (the recombination time 
t rcc ocT ' 7 ). 

We show in Figure^the optical depth to Thomson scat- 
tering r e between the present day and redshift z in these 
four models. These calculations assume that xm = xud 
until z = 3, when helium is instantaneously doubly ion- 
ized throughout the universe. The solid curve corresponds to 
/ion = 0.2 without an early phase of reionization. The short- 
and long-dashed curves correspond to our double reioniza- 
tion model, with z t = 14.5 and Zt = 18.4, respectively. The 
dot-dashed curve corresponds to the model with variable 
/ion (z). T he horizontal dotted line shows the 2<r lower limit 
to r e from lKogut et all d2003l) . We note that, in the double 
reionization models, we can always increase T e by placing 
the onset of the initial reionization phase at arbitrarily high 
redshift. For reference, making the total r e = 0.13 would re- 
quire the fully ionized phase to begin by z ~ 18.4 or z ~ 22 
if z t — 14.5 or zt = 18.4, respectively. 
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3 21 CM RADIATION FROM THE 
INTERGALACTIC MEDIUM 

The optical de pth of the ne utral IGM to redshifted hyperfine 
absorption is 
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Here hp is Planck's constant, k is Boltzmann's constant, 
h>o = 1420.4 MHz is the rest-frame hyperfine transition fre- 
quency, Aio = 2.85 x 10 -15 s _1 is the spontaneous emission 
coefficient for the transition, Ts is the spin temperature of 
the IGM (i.e., the excitation temperature of the hyperfine 
transition), Tcmb = 2.73(1 + z) K is the CMB temperature 
at redshift z, and nm is the local neutral hydrogen density. 
In the second equality, we have assumed sufficiently high 
redshifts such that H(z) ~ Hoftm 2 (l + z) 3 ^ 2 (which is well- 
satisfied in the era we study, z > 6). We have set S equal to 
the local overdensity and £hi equal to the neutral fraction. 

In the absence of collisions and Lya photon pumping, 
the HI spin temperature equals the CMB temperature. In 
this case the IGM will not be visible through the 21 cm tran- 
sition. However, either of these mechanis ms can couple Ts to 
Tk, the kinetic te mperature of the gas llWouthuvsenlll95^: 
lFielJl958tll959bl) . While collisions are only e ffective at ex- 
treme ly high redshifts or in collapsed regions, iMadau et alJ 
lll997l) showed that even a relatively small Lya photon field 
is sufficient to couple the spin and kinetic temperatures. 
ICiardi fc Madaul i2003Tl show that (in their simulations) this 
requirement will be easily satisfied for the redshifts of in- 
terest. Because our simulation has a roughly similar total 
emissivity to theirs, we assume that Ts and Tk are closely 
coupled throughout. 

The visibility of the IGM then depends on its ki- 
netic temperature. If Tk <C Tcmb, the IGM will ap- 
pear in absorption against the CMB. However, the IGM 
will likely be heated shortly after the first sources appear 
through X-ray, photoionizatio n, and shock heating (e.g., 
IChen fc Miralda-Escude1l2003l) . Once heating has occurred, 
the IGM appears in emission at a level nearly independent 
of Ts- For simplicity, and because our simulation does not 
include all of the relevant heating mechanisms (such as X- 
rays), we assume Ts 3> Tcmb throughout. At extremely high 
redshifts, when the first sources are just turning on, the IGM 
may instead appear in absorption, but this does not affect 
the main points of our paper. In the limit Ts S> Tcmb, the 
ob served brightness te mperature excess relative to the CMB 
is llMadau et alJll997fl 
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To apply equation (gj to our simulation, we first choose 



a bandwidth Av and beamsize AO. For reference, these cor- 
respond to comoving distances 



L w 1.2 
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over the relevant redshift range. We then choose a slice of 
depth L within the box and divide it into square pixels of 
width 2R. We randomly select one of the axes of the box 
to be the line of sight and randomly select a starting point 
within the box in that direction. Next, we calculate the to- 
tal mass of neutral hydrogen contained in each pixel. The 
radiative transfer code smooths each particle over the SPH 
kernel and adds the appropriate mass to each cell of the ra- 
diative transfer grid. We multiply the total hydrogen mass 
in each of these cells by the corresponding neutral fraction 
computed in the simulation and then assign this mass to 
the appropriate pixel on the sky. Finally, we divide by the 
mean IGM mass at that redshift and apply equation Q in 
order to find the brightness temperature of the pixel. When 
presenting our results, we always show the average of five 
such slices: because of large scale structure, a single slice is 
not necessarily representative of the box as a whole. Note 
that this procedure neglects peculiar velocities of the gas 
particles, which would shift some particles into or out of the 
frequency range we consider. We discuss the effects of this 
simplification below. 

Although an individual IGM feature is too weak to ob- 
serve, the statistical fluctuations across the sky or in fre- 
quency space may be observable, and they contain informa- 
tion about structure formation and the ionization fraction. 
Equation J3) shows that fluctuations in either the density 
S or neutral fraction xhi of the IGM lead to fluctuations 
in the brightness temperature observed against the CMB. 
While the full simulations are necessary in order to com- 
pute xm, we can estimate the level of fluctuations due to 
variations in the overdensity using linear theory. We wish 
to compute the mass fluctuations inside a cylindrical win- 
dow, with the ax is of the cylinde r along the line of sight 
jTozzi et al.ll200d) . The result is 4 dlliev et al]|2002f) 
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wher e P(k) is the power spe ctrum at the present day (taken 
from lEisenstein fc Hulll998T) . Ji(y) is the Bessel function of 
order unity, D(z) is the growth factor at redshift z (nor- 
malised to unity at the present day), and s(k,x) corrects 

4 We note that our numerical evaluation of equa tion 171 does not 
reproduce the results shown in iTozzi et aflteOOCD : the predictions 
in their Figure 1 exceed ours by ~ 20%, given identical cosmolog- 
ical parameters and an identical redshift space correction. I. Iliev 
has also evaluated the integral and found results similar to ours 
(I. Iliev 2003, private communication). The quantitative match 
with the simulation results shown in Figure|5]gives us confidence 
in our result. 



High-Redshift 21 cm Emission 5 




A0 (arcmin) 



Figure 2. The brightness temperature fluctuations in a neu- 
tral IGM as a function of scale at various redshifts, assuming a 
bandwidth Av = 0.1 MHz. The solid curves show the root mean 
square fluctuation from the simulation assuming a fully neutral 
IGM. The dashed curve shows the estimate from equation 171 
with s ( k, x) = 1. The dotted curve shows the estimate with the 
iKaiserl Il987l> redshift space correction applied to all scales (see 
text) . 

for redshift space distortions. In the limit in whic h L 3> R, 
s(k,x) = (1 + fx 2 ) 2 , where j « fi^ 6 l)Kaiserlll987t ). We note 
here that the IKaiserl 111987ft correction cannot be applied 
outside of this limit. If R> L, the dominant contribution to 
(8Tfr )est comes from modes with ~ R > L, while Kaiser's 
derivation only applies to modes with fc _1 <C L. For R^> L, 
s(k,x) approaches unity. Qualitatively, this occurs because 
the small-fc modes sampled by such a volume are transverse 
to the line of sight and hence are unaffected by peculiar 
velocities. 

We show how our simulation results compare to this es- 
timate in Figure |2] In each panel, the solid line shows the 
simulation results for the root mean square fluctuation as- 
suming a fully neutral IGM, while the dashed line shows the 
estimate of equation J7J with s(k,x) = 1 (i.e., ignoring pe- 
culiar velocity information, as we do in our simulation anal- 
ysis). We see that equation J7J provides a remarkably good 
estimate over this entire redshift range. Note that the sim- 
ulation results steepen slightly with cosmic time, while the 
shape of the estimate is fixed [the only redshift dependence 
is in the growth factor D(z), which is independent of scale]. 
The steepening thus reflects nonlinear structure formation. 
The simulation results begin to drop steeply for A8 > 1', 
indicating that the simulation cannot fully capture fluctua- 
tions on these scales. This is not surprising considering that 
the box subtends only #box ~ 5'. Comparison to equation 
J7J also shows that we underestimate the fluctuations for 
bandwidths Ai/>0.5 MHz. We focus on scales A#<1' and 
Av<0.1 MHz in most of the following. 

Figure [5] also illustrates how the inclusion of peculiar 
velocities would affect our results. The dotted line shows 



the estimate of equ ation l|71 with the linear theory redshift- 
space correction of iKaisen dl987l) applied blindly to all k- 
modes. The fluctuation signal increases by ~ 25%. Of course, 
with our choice of Av, this correction rigorously only applies 
for A8 -C 1', and redshift space distortions will in reality 
vanish as A8 increases. The dotted lines therefore actually 
overestimate the effect of peculiar velocities. We see that 
peculiar velocities will not affect our qualitative results, so 
we chose not to complicate the analysis by including them. 
As a result, we underestimate the fluctuations by < 25% for 
angular scales small compared to the bandwidth. 



4 RESULTS 

Figure[3]shows maps of the IGM brightness temperature 8Tb 
in 21 cm emission at different redshifts. We show the evo- 
lution of a single slice of width 10/i -1 comoving Mpc with 
Av = 0.1 MHz in the f ion = 0.2 model. (Note that the angu- 
lar size subtended by the box varies between redshifts.) The 
nine slices are spaced equally in cosmic time and were cho- 
sen to include the overlap epoch. The colourscale shows 8Tb 
on a logarithmic scale. We see the first HII regions appear 
and grow, eventually filling all of space at z ~ 8. Before 
overlap, fluctuations arise primarily because of variations 
in the density of the IGM. After overlap, filaments stand 
out more clearly than their overdensities would naively im- 
ply. This amplification occurs because recombination occurs 
more quickly in dense regions, so the residual neutral frac- 
tion in filaments exceeds that of voids. During overlap the 
expanding HII regions produce the strongest fluctuations. 

Figure |21 suggests that the mean 8Tb declines dramat- 
ically during overlap. The absolute 8Tb at a single redshift 
cannot be measured because foreground emission (primarily 
the CMB and Galactic synchr otron radiation) gr eatly ex- 
ceeds the IGM signal. However. IShaver et alJ il999Tl pointed 
out that this transition can be observed as a "step" in 
frequency space, provided that overlap occurs relatively 
quickly, because the foregrounds are smooth power-laws in 
frequency. We show in Figure 21 the mean 8Tb as a func- 
tion of z for each of our models. Here the solid line assumes 
/ion = 0.2 and the dot-dashed line assumes a variable ion- 
izing efficiency fi on (z). The short- and long-dashed lines as- 
sume /i on = 0.2 with zt = 14.5 and 18.4, respectively. The 
dotted line shows the mean signal assuming a fully neutral 
IGM. The decline at overlap (z ~ 10-8 in the / ion = 0.2 
models) from 8Tb ~ 25 mK to near zero is obvious. A up- 
ward step caused by recombination after an early ioniza- 
tion epoch will also be clearly visible and of similar mag- 
nitude. Note that the Zt = 14.5 model never fully recom- 
bines, so its 8Tb remains somewhat smaller than the case 
with no early reionization epoch. The decline at overlap is 
less rapid for the variable ionizing efficiency model, both be- 
cause of the shorter recombination time during reionization 
and because (by construction) /i on decreases with cosmic 
time. A given redshift corresponds to an observer frequency 
v = 142 [10/(1 + z)\ MHz. In both the early and late reion- 
ization scenarios, the fluctuations decrease by an order of 
magnitude over a frequency range Av ~ 25 MHz. 

Unfortunately, detecting individual features in Figure 
[3] (rather than the global transition from a neutral to an 
ionized IGM) would require exceedingly long observations, 
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Figure 3. The brightness temperature from the 21 cm transition at various redshifts in the /; on = 0.2 model and with a bandwidth 
Au = 0.1 MHz. Each panel corresponds to the same slice of the simulation (of width 10h~ 1 comoving Mpc), at z = 12.1, 11.2, 10.4, 9.8, 
9.2, 8.7, 8.3, 7.9, and 7.6, from top left to bottom right. 



even with the next generation of radio telescopes. Instead we 
may hope to measure the statistical brightness temperature 
fluctuations caused by variations in the neutral hydrogen 
density. We first consider fluctuations across the plane of 
the sky at a fixed frequency. Figure shows the root mean 
square variation (5T 2 ) 1 / 2 as a function of the beamsize at 
a series of redshifts. In each panel, the dotted line assumes 
a fully neutral IGM, the solid line assumes f lon = 0.2, the 
dot-dashed line assumes a variable fi on (z), and the short- 
dashed line assumes zt = 14.5. All assume a bandwidth 
IS.v = 0.1 MHz. Note that the finite box size artificially sup- 
presses the fluctuations on scales A9 > 1' (see Figure |5J. 

Comparing the four panels, we see that the fluctuations 
remain approximately constant with time until overlap is 
complete (at z ~ 12 for the variable ionizing efficiency case 
and z ~ 8 for the other two cases). Until this time, the fluc- 
tuation amplitude in the models including ionizing photons 
is approximately the same as that of a fully neutral medium. 
At overlap, (ST 2 ) 1 / 2 declines sharply. In the double reion- 
ization scenario the fluctuations are strongly suppressed at 
z > zt but quickly return to near the level of a neutral IGM. 



This is simply because the recombination time is short at 
high redshifts. 

We therefore find that, in principle, measuring the 
magnitude of the fluctuations as a function of scale at a 
fixed frequency is an effective discriminant between dif- 
ferent reionization models. Unfortunately, fluctuations in 
the foregrounds present substantial problems for actually 
measuring IGM fluctuations. While Galactic foregrounds 
will be smoot h on these scales , extragalactic foregrounds 
may not be. iDi Matteo et all (l2002f) extrapolated the 
known population of low-frequency radio point sources to 
smaller flux l evels and estimated {ST^} 1 / 2 > 100-1000 mK. 
lOh fc Mackl |20031) calculated the fluctuations from free- 
free emission by ionized haloes at high redshifts and found 
(STg) 1 / 2 > 100 mK. Either of these signals would be enough 
to swamp the variations we predict before and during over- 
lap, so seeking fluctuations across the sky is in reality un- 
likely to provide much information about reionization. 

Fortunately, fluctuations should be easier to extract in 
pencil- beam surveys of a fixed beamsize made over a large 
frequency range. Here the sampled volume is a long cylinder 
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Figure 4. The mean brightness temperature enhancement as 
a function of redshift for the different reionization models. The 
dotted curve assumes a fully neutral IGM. The dot-dashed curve 
assumes a variable fi on (z); all other curves assume /; on = 0.2. 
The short-dashed and long-dashed curves assume the IGM is 99% 
ionized until zt = 14.5 and zt = 18.4, respectively. 



Figure 5. The mean brightness temperature fluctuations as 
a function of scale at various redshifts, assuming a bandwidth 
Av = 0.1 MHz. Each curve corresponds to a different reioniza- 
tion scenario: a fully neutral IGM (dotted lines), /; on = 0.2 (solid 
lines), a variable /i on (z) (dot-dashed lines), and an ionized IGM 
above zt = 14.5 (short-dashed lines). 



divided into segments by the channel width of the radio re- 
ceiver. Nearly all known foreground contaminants produce 
either free-free or synchrotron radio emission and thus have 
smooth power-law spectra, so fluctuations in the observed 
signal with frequency would be due either to I GM structure 
or to v ariations in the beamsize with frequency. lOh fc Mackl 
( 2003) estimate that the latter effect will (at worst) pro- 
duce fluctuations similar in magnitude to those of the IGM. 
Careful beamsize control will reduce the noise further. 

We therefore show in Figure E] the root mean square 
fluctuations as a function of redshift at a fixed angular scale 
A6> = 0.1' and bandwidth Av = 0.1 MHz. The lines describe 
the same reionization models as in Figure 0] In general, the 
behavior of (ST 2 ) 1 / 2 is similar to that of <5Tj,: the fluctuations 
are large when the medium is primarily neutral and near zero 
otherwise. The transition from a neutral to an ionized IGM 
(or vice versa, in the double reionization models) is thus 
clearly visible as a "step" in the mean level of fluctuations. 

However, there are some subtle but interesting differ- 
ences between the mean signal and its fluctuations. First, 
Figure 2] shows that in a fully neutral IGM 8Tb decreases 
with cosmic time because of the expansion of the universe. 
In contrast, (ST 2 ) 1 / 2 increases with time because of the con- 
tinuing collapse of baryons into sheets, filaments, and bound 
structures. Second, (ST 2 ) 1 / 2 shows a sharper decline than 
the mean signal. The fluctuations remain large during the 



5 An exception would occur if the line of sight intersects a 
galactic disk, which could produce radio recombination line fea- 
tures. However, the se disks should produce very small fluctuations 
iOh fc MaclJl2003l) . and in any case recombination lines have a 
predictable structure. 



early stages of overlap and fall abruptly only at the tail end 
of reionization. We find that, in all scenarios, the fluctu- 
ation amplitude decreases by an order of magnitude over 
Aj/<15 MHz. Intuitively, this occurs because the contrast 
between neutral and ionized regions dominates the fluctua- 
tions during early overlap; until neutral regions become quite 
rare, the fluctuations therefore remain large. In the model 
with a variable /i on (z), (ST 2 ) 1 / 2 exceeds the fluctuations in 
a neutral medium during the early stages of overlap precisely 
because of the contrast between fully neutral and ionized re- 
gions. Finally, note the "noise" in the estimated fluctuations 
after overlap, particularly in the variable /i on (z) model. This 
indicates that, at these times, {ST 2 ) 1 / 2 is more susceptible 
to large scale structure variations beyond overlap (i.e., av- 
eraging five slices is not a sufficiently representative volume 
past overlap). As in Figure|3] this occurs because the shorter 
recombination time in dense regions exaggerates large scale 
structure. Because the post-overlap fluctuations are unob- 
servable even in the best cases, we do not worry about this 
further. 

Figure |7| shows how the signal varies with the choice of 
A6 and Av. Because a larger beamsize or bandwidth implies 
that each pixel samples a larger volume, the fluctuations 
decrease slowly as these parameters increase. However, our 
qualitative conclusions are unaffected by the choice of scale: 
the level of fluctuations remains approximately constant un- 
til the late stages of overlap, when they decline suddenly. In 
fact the decline with scale is not quite as severe as implied by 
this Figure, because the finite box size causes us to under- 
estimate the fluctuations on large scales by a small amount 
(see Figure |5J. 
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Figure 6. The brightness temperature fluctuations as a func- 
tion of redshift in the different reionization models, for a fixed 
angular scale A9 = 0.1' and bandwidth Au = 0.1 MHz. The 
dotted curve assumes a fully neutral IGM. The dot-dashed curve 
assumes a variable fi on (z); all other curves assume /i on = 0.2. 
The short-dashed and long-dashed curves assume that the IGM 
is 99% ionized until zt = 14.5 and zt = 18.4, respectively. 
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Figure 7. The brightness temperature fluctuations as a function 
of redshift for the /; on = 0.2 model. In the top panel, we show 
results for Au = 0.1 MHz and A6 = 0.1', 0.5', 1.0', and 2.0', from 
top to bottom. In the bottom panel, we show results for A0 = 1.0' 
and Au = 0.1 and 0.5 MHz, from top to bottom. In each panel, 
the dotted curves assume a fully neutral IGM while the solid 
curves use the ionization fractions from the simulation. 



5 DISCUSSION 

We have used a high resolution hydrodynamic cosmologi- 
cal simulation with the radiative transfer of ionizing pho- 
tons added in a post-processing step to examine the 21 
cm emission signal from neutral hydrogen at high red- 
shifts. We assume: (1) the existence of a sufficiently strong 
Lycv background to couple the HI spin te mperature Tg 
to th e kinetic temperature T k of the IGM ijMadau et alJ 
Il997t ICiardi fe Madaul 12003ft and (2) a "hot" IGM sat- 
isfying Tk S> Tomb , with heating from, e.g., X-rays 
ijChen fc Miralda-EscudJl2003ft . If these assumptions hold, 
we find a mean excess brightness temperature STb ~ 20- 
30 mK with fluctuations {8Tj*) 1/2 ~ 5-10 mK on small 
scales before reionization. During this era, the mean signal 
decreases slowly with cosmic time because of the expansion 
of the universe, but the fluctuations remain approximately 
constant. At overlap, both the mean signal and the fluctua- 
tions drop rapidly to near zero. Interestingly, while the mean 
signal disappears over a redshift interval Az > 2, the fluctu- 
ations decline only in the final phase of overlap, dropping 
to zero over an interval Az<l. The duration over which 
these two signals fall therefore contains information about 
the process of overlap. 

Combining observations of high-redshift quasars 
jBecker et alJ l200ll : iFan et all 12002ft . which suggest that 
reionization ended at z r ~ 7, with re sults from WMAP 

l Eool . 



ijSoergel et all 120031 iKogut et al J 12003ft . which apparently 
require z r > 14, paints an intriguing picture of reionization. 
Additio nal constraints from t he Lyq forest suggest ing that 
z r < 10 jTheuns et alJ 120021: iHui fc Haimanl 12003ft require 
an even more complex reionization history (although 
these constraints are subje ct to uncertain t ies ab out the 
effects of Hell reionization). ISokasian et al.l <l2003ft showed 
that we can reconcile the quasar and CMB observations 
by smoothly increasing the efficiency of ionizing photon 
production without changing the sources responsible for 
reionization. In such a model, overlap occurs at z> 12, 
but the residual neutral fraction remains large until z ~ 6 
because of the short recombination times. Another way 
to reconcile these observations is wi th two ep o chs o f 
reionization JWvithe fc Loebl l2003allbl ICenl I2002L l2003ft . 
In these models, feedback from the sources responsible 
for the first reionization event halts or modifies their 
subsequent formation and growth. The concomitant sharp 
drop in the efficiency of ionizing photon production al- 
lows some fraction of the IGM to recombine and remain 
mostly neutral until star-forming galaxies at z < 8 ionize 
it again. Because our simulations cannot resolve the small 
haloes responsible for the first era of reionization in these 
models, we have approximated "double reionization" by 
simply fixing xm = 0.01 at high redshift. In the future, 
simulations of reionization by Population III sources will 
allow such scenarios to be studied in more detail (Sokasian, 
Yoshida, Abel, Hernquist, & Springel 2003, in preparation). 
A different class of semi-analytic models, which handle 
feedback-induced self-regulation of ionizing sources in a 
perhaps more sophisticated manner, pre dict a long era of 
parti al ionization between 6 < z < 15 ijHaiman fc Hold"e3 

Although these reionization histories yield identical op- 
tical depths to electron scattering and identical transmission 



statistics in quasar spectra, we have shown that they can be 
easily distinguished through 21 cm observations. An early 
overlap model causes both the mean signal 5Tb and the fluc- 
tuations (ST b 2 ) 1//2 to drop to near zero at z > 13. In a double 
reionization model, on the other hand, both the mean and 
fluctuations return to levels close to that of a fully neutral 
medium after the gas has recombined (see Figures 21 an d 
Ej. Our simulations do not allow us to study self-consistent 
models with long eras of partial reionization. Nevertheless, 
our results suggest a clear signature for such a reionization 
history. The mean signal 5Tb is simply proportional to the 
ionized fraction, so it wo uld decline slowly w i th co smic time 
in a model like those of lHaiman fc Holder! <l2003l) . On the 
other hand, we have found that (5Tb 2 ) 1 ' 2 remains large until 
the tail end of overlap. We therefore suggest that these mod- 
els can be identified by a slowly decreasing mean signal ac- 
companied by large fluctuations. Thus, all of the reionization 
scenarios proposed to date can be distinguished in a straight- 
forward manner through 21 cm emission measurements at 
v ~ 100-200 MHz, corresponding toz~ 14-6. This should 
be contrasted with predictions for advanced CMB polariza- 
tion measurements, in which the d ifferences between these 
sce narios are quite subtle (s ee, e.g., iHolder et aJll2003l) . 

ICiardi fc Madaul i2003T) have recently st udied 21 cm 
emissi on with the reionization simulations of ICiardi et alJ 
(2003). They examined early and late reionization scenarios 
and, like us, found an abrupt drop in (5Tb) 1 ' 2 at the end of 
overlap in both cases. We note, however, that the fluctua- 
tion level they predict exceeds ours (and our semi-analytic 
estimate) by about a factor of two. The most likely explana- 
tion is the differing resolution of our simulations. Our bet- 
ter mass resolution allows us to follow small-scale density 
fluctuations and to identify smaller ionizing sources, which 
signif icantly alter the reionization process iSokasian et alJ 
2003). Moreover, our simulation follows the gas hydrody- 
namics and more accurately d e scribe s gas clumping. On the 
other hand, ICiardi fc Madaul i2003T) have a box twice the 
width of ours, so they are better able to study fluctuations 
on large scales. In any case, this discrepancy does not affect 
any of the qualitative conclusions discussed above. 

Minihaloes in the IGM will also emit 21 cm radiation 
Jlliev et al.ll2002h . Such minihaloes are too small to be re- 
solved by our simulation; however, comparison of our results 
to the semi-analytic model of Iliev et al. shows that emission 
from the IGM will dominate that from minihaloes provided 
that o ur assumption of Tg ~ Tk 2S> Tcmb holds. As pointed 
out bv lOh fc Mackl i2003f) . this is a simple consequence of 
the fact that only a small fraction of the IGM collapses 
into minihaloes. However, the high overdensities in these ob- 
jects (5 > 100) means that the hydrogen spin temperature in 
minihaloes exceeds the CMB temperature even w ithout cou- 
pling from a diffuse Lya radiation background lllliev et alJ 

2002) . Minihalo emission will therefore dominate before such 
a background has been generated (i.e., around the time when 
the very first sources appear). The Lya background will 
rapid ly build up to a point where Ts ~ Tk iCiardi fc Madaul 

2003) , after which emission (or absorption) from the diffuse 
IGM will dominate. During this later era, the most promis- 
ing way to study minihaloes i s through absorption spe ctra 
of high-redshift radio sources jFurlanetto fc LoeblEooal . 

There are two obstacles to observing the 21 cm emission 
signal. The most basic is instrumental noise. The root mean 
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square noise level in a single channel of a radio telescope is 

.arca. . 4 ,„k(!^)'(=^) 

/ 2' \ 2 /0.5 MHz 10 days\ 1 I 2 

{m) {-A^ h —) ' (8) 

where A e s is the effective area of the telescope, T sys is its 
system temperature, Ai/ C h is the channel width, t is the inte- 
gration time, and where we have assumed two orthogonal po- 
larizations. We have used the expected instrumental param- 
eters for the Square Kilometer Array. 6 The Low Frequency 
Array 7 is expected to have A e g/T syB ~ 4.5 x 10 2 m 2 K _1 
in the relevant frequency range. We therefore see that, even 
for these advanced instruments, detecting emission from the 
high-redshift IGM requires large bandwidths and/or large 
beamsizes. Although our simulation box is not large enough 
to predict (5T 2 ) 1 ' 2 on the relevant scales quantitatively, we 
have used a semi-analytic estimate to show that the fluctu- 
ations from density variations decline relatively slowly with 
scale (see Figure We have further shown that this esti- 
mate is quantitatively accurate (to < 30%) over the entire 
range of redshifts and scales we study. Our qualitative con- 
clusions will therefore likely not change on the larger scales 
that can be realistically probed. 

The currently unknown foregrounds constitute the sec- 
ond obstacle to measuring this signal. The mean brightness 
temperature will be swamped by the CMB, extragalactic 
sources in the beam, and Galactic synchrotron radiation. 
Howeve r, because all these sources have smooth power-law 
spectra. [Shaver et al.l (1 199911 have shown that the abrupt de- 
cline in the signal at reionization can be observed with even 
a moderately sized radio telescope. Because we seek only the 
mean signal in this case, the beamsize can be made arbitrar- 
ily large (although a small bandwidth is still desirable for 
learning about the duration of overlap). Assuming that our 
volume is representative, our predictions for the change in 
the mean signal 5Tb will not suffer from errors due to the 
finit e box size. 

iDi Matteo et alJ J2002I) and lOh fc Mackl (12003ft have ar- 
gued that the fluctuations on the plane of the sky caused by 
unresolved extragalactic point sources and free-free emission 
from ionized haloes will exceed those from 21 cm emission 
by at least an order of magnitude. Measuring (ST 2 ) 1 ' 2 as a 
function of angular scale at a single frequency is therefore 
unlikely to be feasible. However, again because all relevant 
foreground sources have smooth spectra, detecting fluctua- 
tions as a function of frequency along a given line of sight 
will be possible, so long as the variation in beam size with 
frequency can be controlled (see lOh fc M ack 2003]). Because 
the fluctuations are relatively constant for transverse scales 
smaller than the length corresponding to the specified band- 
width, the optimal beamsize will have 2R ~ L (see equations 
and El). 

Each of these techniques offers a route to learn about 
reionization that is complementary to quasar absorption 
spectra and CMB polarization observations. Measurements 
of either the mean signal or fluctuations about the mean 



u bee, e.g., http:// www.usska.org/main.html 
7 See |http://www.Tota r.org/index.h tml| 
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will tell us about the history of neutral gas in the IGM, and 
combining the two will tell us about the overlap process. Al- 
though the observations will be technically challenging, the 
reward - the opportunity to distinguish otherwise nearly de- 
generate reionization histories - is great. 
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